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We present the most general parameterization of models of dark energy in the form of a scalar 
field which is explicitly coupled to dark matter. We follow and extend the Parameterized Post- 
Friedmannian approach, previously applied to modified gravity theories, in order to include inter¬ 
acting dark energy. We demonstrate its use through a number of worked examples and show how 
the initially large parameter space of free functions can be significantly reduced and constrained to 
include only a few non-zero coefficients. This paves the way for a model-independent approach to 
classify and test interacting dark energy theories. 


I. INTRODUCTION 

In recent years, cosmological data from experiments 
with exquisite precision (Cosmic Microwave Background 
measurements (Q-i), la supernovae Q , Baryon Acous¬ 
tic Oscillation surveys i) suggest that ^ 96% of the 
matter/energy content of our Universe is in the form 
of an exotic dark sector. Approximately a quarter of 
the dark sector is believed to be weakly interacting Cold 
Dark Matter, while roughly 70% is in the form of a Dark 
Energy component, a substance with negative pressure 
responsible for the current accelerated expansion of the 
Universe. 

The best candidate for Dark Energy is the cosmolog¬ 
ical constant A. The concordance model of cosmology, 
ACDM, is currently the best fit to observations, but it 
comes along with fundamental questions and problems. 
One of them is the coincidence problem, which poses the 
question why the energy densities of the dark sector com¬ 
ponents are of the same order today, when their cosmo¬ 
logical evolution is very different. A possible solution to 
the coincidence problem is a coupling between the Dark 
Energy and the Dark Matter. The introduction of an 
appropriate coupling does not violate observational con¬ 
straints, and it can change the background evolution of 
the dark sector components in order to offer a solution 
to the coincidence problem. 

A plethora of such dark coupling models can be found 
in the literature (see, e.g. [sl-l^l. In most of these mod¬ 
els, the choice of the coupling is purely phenomenological. 
In a recent paper (3^ . we made further progress at the 
level of construction of such models by identifying three 
separate classes of models of dark energy in the form of 
a scalar field {(j)) coupled to Cold Dark Matter (CDM). 

After constructing general models of exotic dark en¬ 
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ergy or modified gravity and checking their mathematical 
and physical viability (for example by identifying funda¬ 
mental problems like ghosts or strong coupling issues), 
one is interested in testing them against the available 
data to see if they might offer a viable alternative to 
ACDM. 

Currently, there is a pressing need of fast and efficient 
ways to rule out and constrain the large number of cos¬ 
mological models available — it would be practically im¬ 
possible to go through each and every one of them indi¬ 
vidually. The Parameterized Post-Friedmannian (PPF) 
approach offers such a framework and has been applied 
to modified gravity theories (see [i^ for a recent 

overview). In this work we apply the PPF approach in 
interacting dark energy theories and demonstrate its use 
through a number of worked examples. In Section HIl we 
go through the PPF basic principles and general formal¬ 
ism, extending it to the case of coupled dark matter/dark 
energy. In Section IIIII we first demonstrate how a few of 
the most well-known phenomenological models in the lit¬ 
erature fit in to this formalism and then we proceed to 
apply it to the general classes of models we presented in 
[^ . We conclude in Section El 

II. FORMALISM 
A. Basic Concepts 

We start by writing the gravitational field equations of 
a theory as 

= SttG , (1) 

where G^^ is the Einstein tensor of the metric 
is the stress-energy tensor of the known forms of mat¬ 
ter (baryons, photons, neutrinos, etc) that are part of 
the Standard Model of particle physics, is the 

stress-energy tensor of (Generalized) Dark Matter and 
represents the stress-energy tensor of all the un¬ 
known modifications to the gravitational field equations 
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that generate the effect of Dark Energy. Such modifica¬ 
tions may be purely due to a Dark Energy fluid or per¬ 
haps due to a modification of gravity. It may be shown 
that any kind of modification of gravity can be put in the 
form III) (see for instance 0). Let us also note that 
although we start with a Generalized Dark Matter which 
may have non-zero pressure and non-zero shear (43 |. we 
shall later on specialize to the Cold Dark Matter (CDM) 
case where both of these quantities are zero. 

The Bianchi identities tell us that the Einstein tensor 
is divergenceless: 


where a is the scale factor, r is the conformal time and 
7 y is the spatial metric, assumed to be flat. The sym¬ 
metries of the spacetime impose that the only non-zero 
components are the energy density p = —T^q and 
pressure P such that T ^ 

The generalized Einstein equations (ED for this ansatz 
give 

= SirGa^ {psM + Pgdm + Pde) (5) 

and 


= 0 . 


( 2 ) 


which in turn implies that 

sume that the standard model particles do not explicitly 
couple to the dark sector so that y= 0. 
This assumption is well justified by observations which 
strongly constraints such couplings [i^. Furthermore, 
a coupling of the evolving quintessence field to baryons 
would lead to time varying constants of nature, which 
are tightly constrained, see and references therein. 

This leaves us with ^ j.{de)^\^ ^ ^ 

neither part is assumed to be individually conserved. 
Thus we have that 


Y7 rp{GDM)fi _ j _ _y7 rp{DE)^s 


( 3 ) 


where the coupling current represents the energy and 
momentum exchange between the dark sector compo¬ 
nents. 

In what follows we aim to parametrise the coupling 
current in terms of metric potentials and their deriva¬ 
tives as well as the scalar modes that are part of the 
stress-energy tensors of the two dark sector components. 
We shall do that in such a way so that the resulting field 
equations contain at most two time derivatives, or equiv¬ 
alently, each dark sector component obeys two first order 
linearized field equations on an FRW background result¬ 
ing from |3]). We shall proceed by considering first a 
FRW background spacetime and find the relevant equa¬ 
tions that describe the dark sector and then consider lin¬ 
ear perturbations about this background spacetime and 
see how this affects the parameterization. Background 
variables will be signified with a “bar” (unless no confu¬ 
sion can arise, e.g. the scale factor a is always a back¬ 
ground variable) while typically all perturbed tensors will 
be preceded by a <5. For instance, we may split Ji, into 
J 1/ J 1/ S J u . 


B. FRW Background 


V? -2- = SttG (Psm + Pgdm + Pde) , (6) 

a 

where H = ^ is the conformal Hubble parameter and 
dots denote derivatives with respect to r. 

Turning now to the coupling current Ji,, the symme¬ 
tries of the spacetime impose that the only non-zero com¬ 
ponent is 


Q = Jo (7) 

while Ji = 0. The function Q{t) is the background cou¬ 
pling function which is for our purposes a phenomenolog¬ 
ically free function. Specific models of a coupled dark sec¬ 
tor will in general result to specific choices of Q{t) (see, 
for example, 1,113, which are two models we present and 
parametrise in Section uni. 

The 1 ^ = 0 component of |3D gives the field equations 
for the evolution of a particular component indexed by 
as 


pi + i'Hpi{l + wi) = siQ, (8) 


where we have defined the equation of state parameter 
for each /-component as wj = Pj/pi and the constant sj 
takes the values 


Si = 


1 

0 

-1 


DE 

SM fields. (9) 

GDM 


C. Linear Perturbations 

1. The perturbed variables 

We now turn to linear perturbations about the FRW 
background. We shall consider only scalar modes. The 
spacetime metric takes the form 

ds^ =—0^(1 + 2'S)dP — 2a^y i^didz"^ 

+ (P [(1 -I- ^h)jij + Diji^) dx^dxP (10) 


Consider a FRW background spacetime described by 
a metric 


( 4 ) 


where 'k, h, and e are four functions of time and space 
(four scalar modes) and 

Dij = ViVj - 


= P [—dP + 'Yijdx^dx^) 


( 11 ) 
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is a derivative operator that projects out the longitudinal, 
traceless, spatial part of the perturbation. 

Let us now consider the perturbed variables of the flu¬ 
ids. These are the density contrast 5 = 5p/p, the scalar 
mode of the momentum, 0, such that Ui = dViO, the di¬ 


mensionless pressure perturbation If = SP/p such that 
(5Tj = and the scalar mode of the shear E such that 

the shear tensor is Ey = DyE. Putting it all together, 
the stress-energy tensor components for a fluid are 

r°o = -p(i + 5) (12) 

r° = -{p + p)v,e (13) 

n = (p + P)V*(0-C) (14) 

T) = p{w + IV)5^j + {p + P)D^^^ (15) 


D. Dark Coupling Parameterisation 

The goal of this article is to parametrise both of the two 
perturbation variables q and S as linear combinations of 
all other perturbations, such as, the fluid variables <5, 0, 
n and E for each fluid, as well as the metric variables 'L, 
C, h and v. This means 12 variables in total for each of q 
and S. However, this linear combination is not entirely 
arbitrary, but must obey certain rules regarding gauge 
transformations. As we shall see, this reduces the num¬ 
ber of effective independent variables to 10 for each of q 
and S. Before proceeding to the parameterization, let us 
briefly discuss gauge transformations. 


2. Einstein and fluid equations 

The perturbed Einstein equations © are 
n (h + 2V^C^-6-H'^'l' + 2V^r] = SttGo^ ^ p/J/, (16a) 

I 

2p + 2n'ii = S'nGa^Y^{pi + Pi)ei, (16b) 

/ 

-h - 2'Hh P mit + &{n^ + 2'Hj 

-v^ (277 - 2 ^ + 2 C + me) 

= 247rGa2^p/n/ (16c) 


1. Gauge transformations 


The metric in HIOI) is in a form which is not gauge-fixed. 
In other words the four scalar modes are not invariant un¬ 
der gauge transformations 5g^i, —>■ Lg^n generated 

by a vector held Parameterized as ^ 

for two scalar modes and the gauge transformations 
of the metric and fluid perturbations will involve combi¬ 
nations of and and their Hrst time derivatives. 

Consider hrst the variables q and S. We hnd that they 
transform as 


1 

q^qP- 
a L 


Qir + (q — 'HQ) (19a) 


and 


and 

-I- C -f P (b -I- 2() -I- ?7 - vk 

= SttGu^ J 2 ipi + Pi)^i (led) 

I 


Turning now to the huid equations, they are obtained 
by perturbing ([3]). To this purpose we dehne the two 
scalar mode perturbations q and S by 

q = SJo V,S = SJ,. (17) 


We hnd 

Si = SwiELSi -f (1 -f wi) 


V^9i-h-V\ 


-snui p — [q- QSi ], 
PI 


(18a) 


and 


Oi 


H{l-2,wi)p 


Wl 


1 -I- W/ 


91 


Hr 


Ip Wl 


Si 


-f-V^E/ P'HP — 


S 


Pi [ 1 -I- w/ 


— QOi 


(18b) 


where the index I runs over all species (and once again 
let us recall that sde = ^ = —Sgdm while S/ = 0 for all 
other species). 


S^SP-Q^t, 

a 


(19b) 


respectively. Thus if we write q and S' as a linear combi¬ 
nation of the metric and huid variables, variables which 
involve in their transformation must combine together 
so that does not appear overall in the transformation 
of the entire linear combination. 

Now the huid variables transform only with the gauge 
variable ^t, i-e. as (dropping the obvious I indices) 


<5 ^ 
6 

n ^ 




37^(1 Pw) — s 


Q 




n 


-Ct, 

a 

1 
a 


Q 


w — 377(1 P w)w -I- sw— 
P 


(20a) 
(20b) 
Ct, (20c) 


while E is gauge-invariant, hence, all four of them are 
allowed to appear in the q and S parameterization. 

However, the metric variables involve ^.l in their trans¬ 
formation. This means that the metric variables must 
combine together so that is eliminated alltogether. 
Following [40| we can hnd three linear combinations of 
the metric perturbations and their hrst time derivatives 
which transform only with the gauge variable ^.t- These 
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are U = h — V^v and V = iy + 2( as well as h+2V^C. The 
latter one is not independent but is equal to U + V^V. 
Thus out of the four metric scalar modes, we are left 
with two combinations, namely [/ and V, which trans¬ 
form exclusively with and 4/ which transforms with 
■ Explicitely, the transformations are 

U ^ U+-?{^T, 

a 

V ^ V 

a 

4- ^ ^ + — . 

a 

Since, q and 4' contain in their transformation, we 
must allow a further metric variable combination which 
does so, but which doesn’t have higher than 2 nd time 
derivatives. The only possibility is the variable U. 

To summarise, we expect that each of q and S can be 
written as linear combinations of the four fluid variables 
(for each fluid) plus the four variables U, V, 4/ and U. 


( 21 a) 

( 21 b) 

( 21 c) 


2. Completing the parameterization 

Following the discussion above, we start from the pa¬ 
rameterization 

q = Ci 4 ' -|- C2V -\- AiU + A2U + A^Sde + A^Sgdm 
+A 50 DE + AqOgdM + A-jAee + A^AgdM 
+AqYide + AiqTigdm ( 22 ) 

and 

S = -l- C4V BiU -j- B2U B^See + B/^Sqem 

AB^Ode + BqOgdm + B’jAee + B^ItcDM 
-\-BgY.EE + BiqTiGDM (23) 

Performing the gauge transformations in ()22[) we find two 
constraint equations, namely 

Ci=Q- mA2 (24a) 

and 

2C2 = Q-HQ-6HAi + 6{H^-H)A2-^A3 

Pde 

-P£EALa, -a,-a,- ^Ay - i££^Ai 24 h) 

PGDM Pde Pgdm 

Likewise, performing the gauge transformations in (1231) 
we find two further constraint equations, namely 

Cs = -&B2'H (24 c ) 

and 

2Ci = Q-6nBi+6{n'^-'H)B2-^B3 

PDE 

_Pgda£^^ _ ^5 _ Bg - ^^B-J - ^^Bs (24d) 
PGDM PDE PGDM 


TABLE I: Gauge invariant variables. 


^ = -\V - \HV 

f = -if/ + u<t! + \{H - n^)v = ^ + 

5^5-^V 

e^e-\v 

p 


The two constraints (I24ap and (I24bl) are then used to 
eliminate Ci and C 2 from ( 1 ^^ while the two constraints 
(I24cl) and (|24dl) are used to eliminate C 3 and from 
(1231) . The remaining perturbations are written in terms 
of the gauge-invariant variables listed in Table U by com¬ 
bining them with V. The result is 

q = i(Q-7fQ)P + Q^- 6 Ail>- 6 A 2 f 

+A36DE + A4SGDM + A^Ode + AqOgdm 
-|-^7nD£; + A^JIgdM + AgT^EE + AiqY^gDM (25) 

and 

S = ]^QV-QB4^-&B2t + B35DE + B45GDM 

-t-B^OuE + BqOgdm + B’jIIde + B^Agdm 
+ Bg'BDE + BiqY^gDM (26) 

Hence, we are left with 20 free functions in total. 


3. Special case: Cold Dark Matter 


From now on we will assume that the Dark Matter 
fluid is completely cold. This automatically means that 
wgdm = Agdm = ^GDM = 0. We shall further make 
the assumption that the Dark Energy fluid has no shear, 
i.e. Tiee = 0. Furthermore, since there is no possibility 
of confusion we shall set wee = w- 

In general, the pressure perturbation Bee would be 
an independent dynamical degree of freedom (see [d^ for 
an explicit model). However, there are many instances 
where Bee is expressed in terms oi 6 ee and 6 ee via 
equations of state such as the Generalized Dark Matter 
model [ij. As in [i^ we shall also assume that the 
pressure perturbation Hd^: is expressed in terms oi 5 ee 
and Oee via equations of state. However, the usual ex¬ 
pression in [ 4 ^ no longer holds, as it does not transform 
correctly under gauge transformations. An expression 
which does is 


Bee = c^See + {cl — cl) 
+p{0c — See) 


3{l + w)n- 


Pde 


Sde 

( 27 ) 
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where and are the (gauge-invariant) effective and 
adiabatic speeds of sound respectively. It may be shown 
that the divergence of the entropy flux is proportional to 
hdb - c 16 de [Hi, hence, the gauge-invariant “relative 
entropy” parameter /i measures entropy transfer to DE 
due to its motion relative to the CDM fluid. The adia¬ 
batic speed of sound is fixed by the equation of state w 
via 


c 


2 

a 


w + 


Q 

PDE 


W 

mii + w)' 


(28) 


Hence, without loss of generality, we may further set A-j 
and to zero. With these choices, the number of free 
functions is reduced to 12. 

We shall further assume the conformal Newtonian 
gauge for which ^ = 0 (so that H = 0). With this 

choice, the gauge-invariant variables we have defined in 
table [J are equal to the conformal Newtonian gauge vari¬ 
ables. 

Let us now re-state the parameterization as well as the 
necessary evolution equations. The two parameters q and 
S are given by 

q = Q^_6Hi$-6H2($-h'H«')-f Hafois 

-\-A 46 c -f A^0ee + Aq9c (29a) 


and 


S = - 6 Bi^- 6 B 2 {^ + n'i>) + B3SDE 

+B 46 C + B^Oee + Bq9c (29b) 

for unknown functions Ai and Bi with i G 1... 6. 

The evolution equations for CDM are 

Sc = V^9c + 3^ + — {QSc-q), (30a) 

Pc 

and 

9^ = -n9c -f 4- -H ^ {Q9c - S) . (30b) 

Pc 

while the evolution equations for DE are 


See = "SwI-LSee + {^ + w) 


V^9EE + 'i‘^ 


— iHIVEE + 


1 


PDE 


[q — QSee] 


(31a) 


and 


9nE = — 


H {1 - 3w) -h 


1 -I- w 


Ode 


Bee 

l-\-w 




1 


PDE [1 + 


— Q9de 


(31b) 


In the following Section we are going to investigate the 
underlying space of models of coupled DM to DE, and 
show how we can construct a “dictionary” of interact¬ 
ing dark energy theories and their PPF correspondences. 
The same method was applied to modified gravity theo¬ 
ries in IH - 


III. WORKED EXAMPLES 

As a “warm-up” exercise, we are first going to demon¬ 
strate the use of our PPF formalism for interacting dark 
energy theories by showing that the functions Ai and Bi 
are severely constrained when one considers specific mod¬ 
els which appear often in the literature. These are the 
“coupled quintessence” model @ , a model where oc 
[HIH elastic scattering of model of Dark Mat¬ 

ter and Dark Energy [H, IHl ■ In table [III one can see 
the list of the models we consider with their coefficients 
displayed. 

Following that, we consider the parameterisation of the 
general classes of coupled theories we constructed in . 
More specifically, in [H| we presented three distinct types 
of models of dark energy in the form of a scalar field 
explicitly coupled to dark matter. We used the pull-back 
formalism for fluids and generalized the standard fluid 
action in order to include a dark coupling. The general 
functional form for the combined dark energy and dark 
matter Lagrangian we considered is 

L = L{n,Y,Z,(j)), (32) 

where n is the fluid number density, Y = and 

As an example, within GR, a quintessence 
field and an uncoupled fluid is described by the La¬ 
grangian L = Y + V {(j)) + f {n). 

We then considered three distinct ways to reduce the 
general function dal]) giving rise to the three Types 
of coupled models which we now want to parametrise. 
These are the Type-1 models where L = F{Y, (j)) + f{n, (j)) 
(the coupled quintessence model @ is a subcase of Type-1 
with the choice F = Y + V{(j)) and / = the Type- 

2 models where L = F{Y,(j)) -|- f{n,Z) and the Type-3 
models where L = F{Y, Z, (j)) + f{n). 


A. Specific models 


1. Coupled Quintessence 


Let us start with the coupled quintessence (CQ) model 
suggested by Amendola @ , which is a specific subcase of 
the Type 1 class of models we presented in [s^. The 
scalar field action for this model is 


S = — 




(33) 


where E(((>) is the quintessence potential. If a constant 
coupling parameter 15 a is assumed, the coupling current 
is found to be (see [H for details) 

Jfi = ( 34 ) 

Writing the scalar field as (p = ^ + (p ioi a, background 
field (j) and perturbation tp, the components of the stress- 
energy tensor for this model are (using expressions from 
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[13) 




-V (35) 

cl- 

1 , 2^^^ 

Q-- 

= -I3aPc4> 

(36) 


3^7f-Q 



for the background, where = 

dV 

dcj) 

, and 


SpDE 



eDE = ^, 

4> 

(37) 

SPde 

= - Vck'P, 


Cs = 1 , P 

= 0 (38) 


for the perturbations. The required coupling parameters 
are found to be 



(39) 

Q%. 

(40) 


4> 

Now we read-off the coefficients. They are: 

Ai = A 2 = Aq = 0 A^ = —— 

1 -I- w 

A 4 = Q A 5 = (41) 

B 5 = Q = 0 

2. Model with oc 


this model has 

Q = 0 (46) 

in the background while 

q = 0 S = (pde + PDE)anDcrD{()c — Sde), (47) 

at the level of the perturbations, with nu the proper 
number density of dark matter particles and ajy the scat¬ 
tering cross section between dark matter and dark en ergy 
(also note that w = const and = 1 in this model) [33} . 
We therefore find that the only non-zero coefficients are 

B5 = —{poE + pDE)anDO'D = —Bq. (48) 

Now we turn our attention to the three general Types of 
models in 


B. Type 1 theory of DM coupled to DE 

Type 1 models are classified in [13 via 

L{n,Y,Z,(j)) = F {¥,(/))+ f{n,(l)). (49) 

For the case where the Dark Matter is CDM we further 
have f{n, (j)) = where a[(j)) is a free function of the 

field 4 >- 

From [13 we have that the coupling current is 


In this model, which was introduced in [13 and [13, 
the energy-momentum transfer vector is parallel to 
the dark matter 4-velocity u^. In our notation we have 
Uf^ = a(l -I- 4', V0c) and 

jfi = rintPc(l + Sc)Ufj,, (42) 

with Fint a local constant interaction rate. The back¬ 
ground coupling function is 

Q = aTintPc (43) 


Afj, - Pct^cj)^ 



q 




Pc^(f)(f)^ 



(50) 


(51) 


while the perturbative coupling parameters q and S are 

q = Q{Sc + dt) and S = Q9c- (44) 

Comparing with our general parameterisation scheme, we 
find that the only non-zero coefficients are 

A 4 = Be = Q (45) 

3. Elastic scattering of Dark Matter and Dark Energy 

This model was introduced in [13 and it considered 
an elastic interaction between dark energy and dark mat¬ 
ter. It is a pure momentum transfer model and its back¬ 
ground cosmology remains unaltered. In our language. 


I 1^2 

where ^ and Note that the expres¬ 

sion for S agrees with the one we recovered previously, 
Eq. (l4n)l . and the expression for q agrees with Eq. ([3^ if 
= 0. This is expected, as the CQ model we studied 
is a sub case of Type 1. 

We now need to express p and ip in terms of fluid-type 
variables. For Type-1 theories, we find from [13 that the 
background energy density and pressure are given by 

PDE = -K Pde = -F (52) 

where we have introduced the function 

K{Y, (j)) = 2 YFy - F (53) 
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which will come in handy below (and also for Tvpe-2). 
Furthermore, the perturbed variables of interest are [39j. 


&DE 


ZKy 

K 





2 _ 2 _j_ 0 \F(I) C^K(^ 

“ ^ ‘^{pDE Pde)'H — Q 


C 


2 

s 


Fy 


^ = 0 


(54) 

(55) 

(56) 


where and similarly for F. The expression for 

c^, above (and since also // = 0), says that if F)/, = 0 then 
the scalar field perturbations are adiabatic. 

To express Lp and Lp in terms of fluid-type variables we 
invert the relations (l56)) to get 


8 Z = - 


K 

ZKy 


8 de 


QjKfj) 


Ode 


ip = (57) 


which are valid for adiabatic and non-adiabatic pertur¬ 
bations. 

Let us now calculate the sought-after coefficients. Us¬ 
ing aSZ = — ip and (1571) we find 



1 -I- w 


See 


(1 -I- w)K 


4 > 6 de- 


(58) 


so that (EU gives the required coefficients as 


Ai — A 2 — Aq — 0 A;} 


A 4 = Q 
B^ = Q 


A 5 = 




Qc 


2 

s 


1 -I- w 


{l + w)K 


(59) 


For the case where = I3a is a constant and furthermore 
F{¥,(/)) = y -|- V{cj)) we recover the coefficients for the 
coupled quintessence model discussed above, which is in 
fact a subcase of a Type-1 model of coupled dark energy. 


variables given by (|56l) are still valid for the case of Type- 
2 theory and the function K is still defined via (1551) . 

In order to proceed further we need the function Q 
which at first glance using dSB) is given by 


Q 


z 


{pc + SUpc) p + PcPzZ 


(63) 


Using ([5]) to eliminate the terms pc and Z we find 


^ ZPz ^ 

where Z is determined from eq (70) of as 
SZFy'H -(- aKrf, 


(64) 




Ky--^ 


(65) 


i+zp 


The perturbative variables q and S are also found from 
den. Firstly S is easily calculated as 


S = Q 6 ^ 


( 66 ) 


while q is found to be from (|6T1) as 


:{(i 


q = + {S,-y^9,-3^) pcP + PcPzSZ 


{p, + 3-Hpc) P + PcPz^ [<5c - 4-] 
{pc + ‘iB.pc) Pz + PcPzzZ F 


5Z\ (67) 


However, using (|63)) as well as (l30l) in order to eliminate 
5c the expression for q simplifies to 


q — QSc + Q—r + -T 7 p 
Z dZ 


ZPz 


i + zp\ 


PcZ 5Z 


( 68 ) 


C. Type 2 theory of DM coupled to DE 


Type 2 models are classified via 

L{n,Y,Z,p) = F{Y,p)+f{n,Z), (60) 


with / = nh{Z) in the case that the scalar field is coupled 
to CDM. The coupling current in this case is 


{PcPu’') Vfj^p 


where P{Z) is the function 


P{Z) 


hz 

h — Zhz 


(61) 


(62) 


and hz = ^ (and the same when Z is used as a subscript 
for P). Let us first note that the relations for the fluid 


What remains is now to eliminate 5Z. This can be done 
using the perturbative version of eq (70) of which 
gives 




= If -h < ZKyy + Pc 


Pz 


dZ \1 + ZP 


Ky - AAz_ 
^ i+zp 


ZKyU + Ky^P 


■SZ- 


3ZF. 


Y 


+z 


3ZFyT~L + aK^ J 3ZFyT~L + aK^ 
iFYfPHp — — FyV^ 


($ + 774') 


SZFyB. + aK(j) 


aK 


Y<t> 


Ky - 


PcPZ 

1+zp 


PcPz 


(1 + ZP)Ky - PcPz 


(69) 


Using dMD and (|571) into (1551) we may now determine the 






































coefficients. They are 


Ai — Aq — 0 
A3 = 


Ay, — 


ZFy 


(1 + w) ( 


Ky - 


6 ZFyT~L + 2 aK^ 

c 


Q 




ZKy^Iu ^ 

dz yi + zp 


+Z Kyy 


i+zp J 

p,Pz [q + {iKYU + Ky^$) 


1 + zp 




Af, — 


(1 + ZP)Ky - PcPz 

xK(pQ 


^ Ky - 

^ i+zp 

d ( Zfdz 

dZ^'^yi + Zjd 


PcPz 


+ 


ZKyy 

ZKy{1 + ZP) Ky 
ZpcPzi^KYH + Ky(I)^) aK^ 


1 + zp 

Z^p,Pz{5FY^n$ - - FyV^) 

i + zp 

= Q Bi^3 = 0 


QZaKYcj) \ 


D. Type 3 theory of DM coupled to DE 

Type 3 models are classified via 

L(n,r,Z,0)=F(y,Z,<^) + /(n). (71) 

The coupling current in this case is [s^ 

J,=q^,i^Xyp(^ + FzypZFZFzu^^V^upY (72) 

where X = V^{Fzu^) A straightforward calculation 
gives 


Q = (? = 0 


(73) 


(although there can be 2nd order corrections to Jq) 
This means that the Type 3 case provides for a pure 
momentum-transfer coupling up-to linear order in per¬ 
turbation theory. 

To proceed to the coefficients we need S which is found 
to be 


S = - {X4> + FzZ + ZFzU) 0, - ZFzOc - -Fzp + Xtp 

(74) 

where the background value of X is 


x = i 

a L 


{ZFzy — Fzz)Z — Fz^pcj) — STLFz 


(75) 


We eliminate the 9c term using (l30l) to get 


where to remind the reader SZ = —^ ^ip — . Now we 

need to express ip and p in terms of the fluid variables. 
From we find 

pDE = Z^Fy — ZFz + F Pde = —F (77) 

for the background variables while 

SpDE = Z [Fy — Z'^Fyy + 2ZFyz — Fzz\ SZ 

+ [Z'^Fy,/, — ZFzcj, + F^] p (78) 

SPde = {ZFy - Fz)SZ - F^p (79) 

+ FzBc 


Ode = 


Fz — ZFy 


(80) 


for the perturbations. For completeness, the adiabatic 
sound speed is 


= 


(70) 


SnjZFY -Fz)-a [F^ + Z^Fy^ - ZFz^] 
2i'HZ{Fy + 2ZFyz — Z'^Fyy — Fzz) 
aF^ 


m{ZFY - Fz) 

while the effective sound speed is 
ZFy — Fz 


cl = 


(81) 


(82) 


Z (^Fy + 2ZFyz — Fzz — Z^Fyy) 

and the relative entropy parameter is 

P = - cl){pDE + Pde)'H (83) 

Clearly if = 0 then the perturbations are adiabatic, 
i.e. and ^ = 0 (so that Ade = c^^de)- 


We can now proceed to find the coefficients. Equations 
(1781) and (l80l) can be inverted to give 


5Z = 


_P_ _ aF^ 
Fz Fy 
ciZ 


Ode + _ 9c 

PDE + E^DE 


1 -I- W 


5de 


and 


, Fz „ aFz 

p = a { j;- - Z j Ode - 


(84) 


(85) 


We also need the equation for Z which is found to be (eq 
(75) in [13 ) 


Pc 


Xp - {X^ + FzF^c + Fz5Z 

■ (76) 

Z - -iUZ 

- 




[“ m{ZFY-Fz)\ 
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Model / Coefficients 

Q 

Ai 

742 

A 3 

A 4 

As 

746 

Hi 

B 2 

Bz 

Ba 

Be 

Be 

Coupled Quintessence 

I^Apc^ 

- 

- 

Q 

1+tt; 

Q 

PaPcO^V^ 

- 

- 

- 

- 

- 

Q 

- 

tJ fj, OC Xlfx 

Q-Pint pc 

- 

- 

- 

Q 

- 

- 

- 

- 

- 

- 

- 

Q 

elastic scattering 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

—Pdb(1 + w)ano(yD 

-Be 

Type-1 

pcO/-(j)(j) 

- 

- 

qA, 

1+tt; 

Q 

<^ 4 . (l + m)K 

- 

- 

- 

- 

- 

Q 

- 

Type-2 

Zpzpc ^ 
1+Z8 ^ 

- 

742 

743 

A 4 

As 

- 

- 

- 

- 

- 

Q 

- 

Type-3 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Bz 

- 

Be 

D , snzFzci 

' Z Fz 

_ Pc _ 


TABLE II: Specific models and their PPF coefficients. The coupled Quintessence model is a subcase of Type 1 with = Pa- 
The elastic scattering model is in fact distinct from Type-3 (see text at the end of section llllDll . For the coefficients A 2 , A 3 , 
A 4 and As in the case of Type-2 see (I70II . For the coefficients Bz and As in the case of Type-3 see (I 86 II . For the remaining 
functions the reader is referred to each specihc example in the text. 


The above equations are then inserted into (TTSl) to give 
the required coefficients as 


Bi 

B 3 

B 5 

Be 


B 2 — B^ — 0 
1 ZFzcl 

1 _ 1 -p ry 


a 

\x(^^ z\ 

+ Fz 


F]] 

2 _ ZFz 

[ Ur J 

aFz 

Fy\ _ 


—Be -\- 


mZFzcl 

_ ZFz 

pc 


( 86 ) 


It would seem tempting to try model the elastic scat¬ 
tering model [13 discussed above (section IIIIA3I) into 
the Type 3 class. However, this is in fact impossible. As 
we can easily check, the elastic scattering model requires 
Be = 0. Within the Type 3 class this is possible only 
if F is independint of Z (i.e. Fz = 0). This implies 
that Be and Be are also zero, in other words, the model 
becomes completely uncoupled. Hence, it is impossible 
construct a model of elastic scattering between CDM and 
DE within the Type-3 class of coupled Dark Energy. 


IV. CONCLUSIONS 

We have presented the most general parameterization 
of models of Dark Energy which is explicitly coupled to 
Dark Matter using the Parameterized Post-Eriedmannian 
framework, and have shown that it is able to encapsulate 
a rich variety of theories. 

Starting from the linearised Einstein equations and us¬ 
ing the Bianchi identities we managed to express the 
modifications to GR coming from the dark sector cou¬ 
pling as a collection of new terms containing the met¬ 
ric potentials and their derivatives as well as the scalar 
modes of the two dark sector components, i.e. the fluid 
variables of (generalised) Dark Matter and Dark Energy. 
Of course, our formalism is based on a few basic as¬ 
sumptions: the background cosmology has an FRW solu¬ 
tion, all field equations are at most second-order in time 


derivatives, and the field equations are gauge-invariant. 
Completing the parameterization we were left with 24 
free functions, but demanding gauge invariance we de¬ 
rived 4 constraint equations which eliminated 4 free func¬ 
tions. 

Twenty free functions in our general parameterization 
is certainly a big number, but by imposing certain well 
motivated assumptions, for instance that the Dark Mat¬ 
ter is Cold, that the Dark Energy is shear-less and that 
the pressure perturbation is not a dynamical quantity we 
reduced the number of free functions to 12. Furthermore, 
we showed that only a handful of these functions are 
non-zero when one considers known models. We demon¬ 
strated this by investigating a number of specific models 
in the literature, as well as the classes of theories we con¬ 
structed in l3M . It is useful to note that, although our 
theories in |39| are derived from an action, the PPF pa¬ 
rameterization does not require knowledge of the action, 
but only knowledge of the field equations. This means 
that the PPF parameterization is a very useful tool for 
phenomenological model building (see [42l | for further dis¬ 
cussion in the context of modified gravity theories). The 
full list of models we consider in this work is displayed in 
table [H] along with their coefficients. 

Our Type 1, 2 and 3 theories contain a fairly gen¬ 
eral coupling function and hence they encapsulate many 
different models. The parameterization coefficients for 
these theories can depend, of course, on the chosen cou¬ 
pling function and its derivatives, and other quantities 
such as the background coupling Q, the background field 
energy density the quintessence potential V{(j)), the 
speed of sound etc. For Type 1 theories there is only 1 
non-zero B coefficient and 3 non-zero A coefficients, for 
Type 2 there is 1 non-zero B coefficient and 4 non-zero 
A coefficients, while for Type 3 all A’s are automatically 
zero and there are 3 non-zero H’s: different classes of 
theories correspond to different non-zero functions. In 
particular, from all the cases we have studied, the coeffi¬ 
cients Ai, Ag, Hi, B 2 and H 4 are always zero. It would 
be indeed very interesting to find models for which any 
of these coefficients is non-zero. 

It would also be interesting to consider the inverse 
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problem, i.e. given Q{t) and a set of PPF coefficients 
Ai and Bi can we reconstruct the functions that appear 
in the coupled Dark Energy Lagrangian for each of the 
Type-1,2 and 3 theories? For instance, what kind of func¬ 
tions correspond to constant PPF coefficients? Tackling 
the inverse problem will help to reduce the free functions 
into simple functional forms which are parametrised by a 
set of constants and can make constraining such theories 
easier and more efficient. 

Another important question is how we can further con¬ 
strain the PPF functions? As discussed in , we might 
expect to find that a subset of the PPF functions can 
be very well constrained, whilst another subset can not. 
However this might not pose a serious problem, as the 
constraining power of a few PPF functions might be suf¬ 
ficient to distinguish between theories. Tackling the in¬ 
verse problem will certainly help here as it can guide us 
to which functional form of the PPF coefficients is the 
most useful. The implementation of the PPF framework 
presented here in numerical codes for the computation of 


the cosmological effects of interacting Dark Energy could 
provide an answer to these questions. We plan to inves¬ 
tigate this in future work. 
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